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ABSTRACT 

Recent observations have allowed the geometry and kinematics of the M87 jet to be 
tightly constrained. We combine these constraints with historical Very Long Baseline 
Interferometry (VLBI) results and the theory of synchrotron self-absorbed radio cores 
in order to investigate the physical properties of the jet. Our results strongly suggest 
the jet to be dominated by an electron-positron (pair) plasma. Although our conserva- 
tive constraints cannot conclusively dismiss an electron-proton plasma, the viability of 
this solution is extremely vulnerable to further tightening of VLBI surface brightness 
limits. The arguments presented, coupled with future high-resolution multi-frequency 
VLBI studies of the jet core, will be able to firmly distinguish these two possibilities. 
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1 INTRODUCTION 

Jets are intimately linked with accretion processes. They 
can be found in almost all astrophysical situations in which 
matter is believed to be undergoing disk accretion onto 
some central collapsed object. Scales range from the slow 
(~ 100 kms" 1 ) sub -parsec jets associated with forming stars 
to those found originating from powerful active galactic nu- 
clei (AGN) which are initially relativistic and can propagate 
for hundreds of kpc before disrupting. Despite intense study, 
many basic questions remain open. The underlying forma- 
tion mechanism is still uncertain (although hydromagnetic 
processes associated with the accretion disk seem a promis- 
ing candidate mechanism: Blandford & Payne 1982; Lynden- 
Bell 1995 and references therein). Moreover, whereas jets 
from forming stars are believed to be molecular material 
entrained by a faster outflow of atomic material, the basic 
nature of AGN jets has not been unambiguously determined. 
Possibilities include an electron-proton plasma (i.e. normal 
material), an electron-positron plasma (i.e. a pair plasma), 
Poynting flux or combinations of these. This is an observa- 
tionally difficult issue to address: even the handful of well 
studied jets only show power-law continuum radiation (pre- 
sumably of synchrotron origin) with no sign of any spectral 
features from material moving with the bulk flow. This is 
to be contrasted with the Galactic source SS 433 in which 
the bulk flow of the jet is a known source of atomic line 
emission. 

Celotti & Fabian (1993) have addressed the issue of the 
matter content of jets in a sample of radio- loud quasars and 
radio galaxies. Combining the synchrotron self-Compton 



constraints with indications of the kinetic luminosity sug- 
gest that, for the sample as a whole, the jets are either 
cold electron-positron flows or electron-proton flows with an 
electron low-energy cutoff of ~ 50MeV. They argue for the 
electron-proton case on the basis of annihilation constraints 
under the assumption that the pairs originate from the in- 
ner regions of the accretion flow. However, the issue of the 
matter content of any one given extragalactic jet (as op- 
posed to a large, possibly heterogeneous, sample) has not 
been addressed in detail. An (observationally based) deter- 
mination of the matter content of an individual well studied 
jet would be an important step in the study of jet formation, 
propagation and emission. 

M87 (NGC 4486) is a giant elliptical galaxy near the 
centre of the Virgo cluster of galaxies. It is associated with 
the radio source (Virgo-A, 3C 274) and classified as a FR-I 
radio galaxy (Fanaroff & Riley 1974) on the basis of its edge- 
darkened morphology and low radio luminosity (Pi78MHz ~ 
1 x 10 32 ergs -1 Hz -1 ). This galaxy contains the most spec- 
tacular example of an extragalactic jet in the northern sky. 
It was the first extragalactic jet to be discovered (Curtis 
1918) and, since then, it has been subjected to intense ob- 
servational studies at all wavelengths (see Biretta 1993 for 
a recent review). The complex knotty structure of the jet 
has been resolved at wavelengths from the radio through to 
X-rays (Biretta, Stern & Harris 1991). The spectrum and 
proper motions (Biretta, Zhou & Owen 1995) of the knots 
yield direct constraints on the physical processes operating 
and the geometry of the system. Hubble Space Telescope 
(HST) studies of the base of the jet reveal a rotating gas 
disk apparently lying normal to the jet direction (Ford et al. 
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1994; Harms et al. 1994). The projected appearance and ve- 
locity of the disk suggests it to be at an inclination of 42 ±5° 
and to enclose a mass of (2.4 ± 0.7) x 10 9 M . Also, Very 
Long Baseline Interferometry (VLBI) has resolved structure 
within the core of M87 on spatial scales of ~ 0.01 pc. The 
VLBI studies will be discussed in more detail in Section 3.1. 
The synchrotron spectrum of this jet has been recently ex- 
amined by Meisenheimer, Roser & Schlotelburg (1996). This 
wealth of data makes M87 the most promising candidate for 
any physical study of extragalactic jets. 

In this paper we combine these observational con- 
straints of the M87 jet with the theory of synchrotron radia- 
tion in order to constrain the physical parameters of the jet. 
In particular, we address the particle density, magnetic flux 
density and the matter content of the jet. In order to keep 
this paper focussed, our treatment is as model independent 
as possible. The arguments we use depend only upon an ob- 
servationally well determined kinematic model and the stan- 
dard synchrotron jet model of Blandford & Konigl (1979). 
Our arguments are independent of any models (of which 
there are a plethora) for the initial formation, collimination 
and acceleration of the jet. Section 2 briefly reviews the rel- 
evant observations and the kinematic model to which they 
lead (Biretta 1993; Biretta, Zhou & Owen 1995). Section 3 
discusses the synchrotron self-absorption model for the ra- 
dio emission from the core and the subsequent constraints on 
the particle density and the magnetic flux density. Section 
4 places further constraints on the particle density of the 
jet using the large scale limits on the total kinetic luminos- 
ity. The resulting constraints suggest the jet to be electron- 
positron dominated rather than electron-proton dominated. 
Section 5 discusses the robustness of the result and some im- 
plications. It is found that an electron-proton jet possessing 
an electron population with a low-energy cutoff ~ 10m e c 2 
also satisfies our constraints although this solution is ex- 
tremely vulnerable to further tightening of the VLBI sur- 
face brightness lower limit. Section 6 summarizes our con- 
clusions. 

Stellar surface brightness fluctuations give a distance 
to M87 of 15.9 ± 0.9 Mpc (Tonry 1991), independently of 
Hubble's constant. In the rest of this paper we assume a 
distance to M87 of 16 Mpc. 



2 KINEMATIC MODEL 

In this work, we use the kinematic model for the M87 jet 
presented in Biretta (1993; hereafter B93) and Biretta, Zhou 
& Owen (1995; hereafter BZO). This is based on detailed 
observations of proper motions within the knotty structure 
of the jet. In this model, the inner jet is relativistic with 
a bulk flow Lorentz factor, T > 3. The prominent knot A 
marks the location of a strong relativistic shock in which 
the flow slows from velocities of ~ c to c/3 (in the frame of 
reference in which the shock is stationary). In the observers 
frame, the Lorentz factor decreases across the shock from 
r ~ 3 to T ~ 1.5. The flow suffers another shock at knot C 
at which it is decelerated to subrelativistic velocities. The 
jet then disrupts and merges into the radio halo. The ori- 
entation of the jet axis to the line of sight, 9, is between 
30-40° . 

The relatively large angle to the line of sight (30-40°) 



puts severe limits on the extent to which relativistic beaming 
can be important. This fact will be extensively used in the 
later sections. The effects of beaming are characterized by 
the beaming parameter, 8, defined by 



r(i - /?cos6>) ' v ; 

where f3 c is the bulk velocity of the material in the observers 
frame. For a given 9, this has a maximum value of 

5max = sin0' (2) 

corresponding to /3 = cos9. Evaluating for 9 = 30° gives 
<5max = 2 (and increasing 9 will result in decreasing <5 ma x)- 
Thus, relativistic effects can only produce a moderate in- 
crease in the observed luminosity of features moving along 
the jet. Moreover, for T > 7 (given 9 — 30°) relativistic aber- 
ration diminishes the observed luminosity since it is beamed 
away from the observer. 

We note that relativistic beaming is important in ren- 
dering the counter-jet undetectable. Assuming the system 
to contain two intrinsically symmetric (back-to-back) jets, 
a bulk Lorentz factor of F > 3 and 9 = 30° results in a 
jet/counter-jet flux ratio of > 300. The observational lim- 
its on the jet/counter-jet (radio) flux ratio (based on the 
non-detection of the counter-jet) are ~ 150 (Biretta 1993). 



3 SYNCHROTRON SELF-ABSORPTION 
CONSTRAINTS 

3.1 Basic model for the core emission 

Like many extragalactic radio sources, the core of M87 has 
a flat radio spectrum (see Fig. 5 of Biretta, Stern & Harris 
1991; hereafter BSH). This is generally believed to be due 
to the superposition of many synchrotron emitting compo- 
nents, each of which has a different self-absorption frequency 
(Blandford & Konigl 1979). In the diverging flow of the in- 
ner jet, these components correspond to emitting regions at 
different distances along the jet. The observed flux density 
at a given frequency is then dominated by the region of the 
jet which becomes self-absorbed at that (Doppler shifted) 
frequency. 

Consider one such region of the jet which is just becom- 
ing self-absorbed at the jet rest-frame frequency v' m . Radia- 
tion with this frequency will suffer relativistic Doppler shifts 
to an observer-frame frequency v m — v' m 5. Suppose the cor- 
responding flux density from the self absorbed region in the 
observers frame is F m . Observations of the core at the fre- 
quency z/ m will measure this flux density plus contributions 
from optically thin parts of the jet. The total observed core 
flux density at v m will be Ftot(^m) = r\ F m where rj > 1 al- 
lows for the emission from the optically thin parts of the 
jet. Within the flat spectrum regime, most of the emission 
must be dominated by the self-absorbed region of the jet: 
thus we impose the condition that r\ < 2. If the observations 
also resolve the flat spectrum core, the resulting dimension 
is a measure of the diameter of the jet at the self-absorption 
point. 

VLBI observations of M87 have been performed at a 
number of frequencies. Reid et al. (1989) present 1.66 GHz 
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VLBI observations of the nuclear region of M87 show- 
ing a well collimated knotty jet extending from an unre- 
solved core. These observations also show subluminal mo- 
tion within these knots suggesting that the pattern speed 
is very much slower than the fluid speed. The 5 GHz obser- 
vations of Pauliny-Toth et al. (1981) just resolve the core 
to have a (major axis) angular diameter of 8d = 0.7 mas 
(corresponding to 0.06pc at the distance of M87) and a to- 
tal core flux of F m = l.OJy. Higher frequency observations 
probe deeper into the core of the jet before encountering 
the self-absorption surface as well as achieving a better an- 
gular resolution. Thus, 22 GHz VLBI observations with a 
resolution of ~ 0.15 mas (corresponding to O.Olpc: Spencer 
& Junor 1986; Junor & Biretta 1995) reveal an unresolved 
core (peak surface brightness 0.34 Jy beam" 1 ) and a station- 
ary pattern of knots stretching for 4 mas in the approximate 
direction of the large scale jet. The best resolution has been 
achieved with 100 GHz VLBI (Baath et al. 1992) which re- 
veals several unresolved knots within 0.1 mas of the core. 
The flat spectrum nature of the core is consistent with it 
being self absorbed up to frequencies of at least 22 GHz and 
possibly even 100 GHz (Spencer, private communication) . 

In the remainder of this section, we use the 5 GHz VLBI 
constraints on the flux and dimensions of the self-absorbed 
region to place constraints on the relativistic (synchrotron 
emitting) proper particle density, n, and magnetic flux den- 
sity, B, at this location within the jet. We assume a uniform 
jet in the sense that n, B and f3 are constant across cross- 
sections of the jet. 



3.2 Surface brightness of the self-absorbed region 

Here we follow the work of Ghisellini et al. (1992; hereafter 
G92). Suppose the relativistic electrons responsible for the 
synchrotron radiation have a differential Lorentz-factor dis- 
tribution of JV(7) = iVo7 _( - 2a+1 - ) for 7 m i n < 7 < 7max where 
7 is the electron Lorentz factor as measured in the rest frame 
of the jet material. 

In the optically-thin regime (and well away from the 
low-frequency or high-frequency cutoff), this produces syn- 
chrotron emission with a spectral index a. For the M87 jet, 
a ~ 0.5. The proper relativistic particle number density is 
given by 



N(-y)dy. 



(3) 



For a > and 7max S> 7min, this integral is dominated by 
the value of 7mm- The particle number density, n, is then 
related to the normalization of the relativistic electron dis- 
tribution, No, by 



No 



7n 



(4) 



where we have used a — 0.5 in order to obtain the last ex- 
pression. In this case the mean Lorentz factor of an electron 
is given by 



(7) = 7n 



,ln 



7n 



7n 



(5) 



For 7min = 1, any reasonable choice of 7 max gives {7) ~ 10. 

The synchrotron flux in the self-absorbed regime is in- 
dependent of particle density. Thus, our constraints on _F m 



and #d can be translated into a constraint on B. In par- 
ticular, there is a relationship between B and the surface 
brightness of the self-absorbed region (e.g. see G92) which 
gives, 



B < W- 6 b(a)e 4 d F~\ 



1+2 



(6) 



where 8d is in milliarcsec, F m is in Jy and b(a) is a function 
tabulated in G92. The inequality results from the limit on 
5. From G92, b{a = 0.5) « 3.2, and noting that c5 max « 2 
gives B < 0.2 G. This constraint is plotted on Fig. 1 (shaded 
areas represent forbidden regions within the n,B plane.) 



3.3 Optical depth of the self-absorbed region 

The point at which the jet is just becoming self-absorbed at 
an observed frequency ^ m is defined by /t(z/ m )JC = 1 where 
k(v) is the absorption coefficient for synchrotron absorption 
at frequency v and X is the path length of the line of sight 
through the jet. Using the standard result for the rest-frame 
absorption coefficient (Jones, O'Dell & Stein 1974; Rybicki 
& Lightman 1979) and noting that v m — v' m 5 gives 



fc(z/ m ) = 



3 a+1 iv 1/2 g(p)e 2 No 



8m c c 



(3/2+a) -(5/2+a) A (5/2+a) 



(7) 



(cgs units) where vs, is the cyclotron frequency, p — 2a + 1 
and g(p) is the given in terms of Gamma functions by 



g(p) = 



r((3p + 22)/12)r((3p + 2)/12)r((p + 6)/4) 

r((p + 8)/4) 



(8) 



This expression for n(y) is valid for v S> fmin ~ 7min^B, 
where i-w n is the low-energy cutoff in the synchrotron spec- 
trum corresponding to the low-energy cutoff in the electron 
energy distribution. In this regime the self-absorption de- 
pends only on the normalization of the relativistic electron 
energy distribution and not on details of the low-energy cut- 
off. The spectrum of BSH shows no evidence for such a low- 
energy cutoff thereby giving us confidence that the observed 
frequencies lie in the domain of applicability of eqn (7). 

We also need to relate the observer-frame line of sight 
path length through the jet, X, with the physical jet diam- 
eter, 2r. Taking into the account the relativistic transfor- 
mations, the required relation (for a cylindrical geometry) 
is 

X = dmT < 9 > 
Evaluating the above condition on the optical depth for 
a — 0.5 gives nB 2 > 27~ i 1 n (5~ ax (cgs). For the moment we 
assume that 7 m i n = 1 (i.e. that the electron energy distribu- 
tion extends down to non- relativistic energies). Relaxation 
of this assumption will be discussed in Section 5.3. The as- 
sumed geometry implies 8 < 2 leading to nB 2 > 0.5 (cgs). 
This is our second constraint on the n,B plane (Fig. 1). 



4 KINETIC LUMINOSITY CONSTRAINTS 

The total kinetic luminosity of the jet is an important ob- 
servational parameter. In this section we estimate the total 
(time-averaged) kinetic power of the jet by examining the 
energetics of the system as a whole. We then assume the 
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Figure 1. Constraints on the B-n plane imposed by synchrotron self-absorption and total kinetic luminosity considerations. The 5 GHz 
VLBI surface brightness constraint on B gives line A. Line B corresponds to the condition that resolved core of the jet is just becoming 
optically thick to self-absorption. The horizontal lines are the densities derived by imposing L43 = 1 for the two cases of a pair jet (with 
7min = 1 and (7} RJ 10) and a cold jet of normal matter. 



VLBI jet carries this same kinetic power in order to place 
independent constraints on its properties. 

4.1 The global energy balance 

Ultimately, the kinetic energy carried by the jet can be ra- 
diated (at the knots, for example), can perform work on the 
surrounding interstellar medium (ISM) to expand the radio 
halo or can be stored in relativistic particles within the ra- 
dio halo. Here we assess the kinetic luminosity of the jet by 
summing these components. 

4-1.1 Radiated energy 

A multiwaveband study of the radiant emissions from the 
jet has been performed by BSH. By summing the observed 
emission of the jet (dominated by emission from knots A and 
B) as tabulated in BSH, we estimate that ~ 3x 10 42 ergs -1 
is radiated along the length of the jet. As stated above, rela- 
tivistic beaming will not have a major effect on the observed 
luminosity. 

4-1.2 Energy stored within the radio halo 

The large scale radio morphology of the system (Turland 
1975) suggests the presence of an inner radio halo (~ 5kpc 



in extent) embedded in a larger, more diffuse, halo (~ 50kpc 
in extent). The inner structure is clearly associated with the 
current phase of nuclear activity since the jet can be obser- 
vationally traced from its origins in the nucleus to the point 
at which it disrupts and feeds this radio halo (e.g. see VLA 
images of Owen, Hardee & Bignell 1980 and, more recently, 
Hines, Owen & Eilek 1989). Spectral curvature arguments 
suggest an age of ~ 10 yr for this inner structure. The ex- 
tensive outer halo is probably associated with a previous 
phase of activity and will not be considered further here. 

Estimating the energy stored within the inner radio halo 
amounts to determining the pressure of this structure. Radio 
observations allow the minimum (i.e. equipartition) pressure 
Pmin of the synchrotron emitting plasma to be determined. 
The VLA data presented by Hines, Owen & Eilek (1989) 
show the inner radio lobe to have a filamentary structure 
with the minimum pressure of the average filament being 
Pmin ~ 2-6 x 10- 10 dyn cm . The brightest filaments have 
Pmin ~ 10~ 9 dyncm~ 2 whereas the interfilamentary plasma 
has Pmin ~4x 10 -11 dyncm -2 . X-ray observations of ther- 
mal emission from the interstellar medium (ISM) also allow 
the ISM pressure to be determined at the vicinity where it 
comes into contact with the inner radio halo. A deprojec- 
tion analysis (Fabian et al. 1981) of data from the ROSAT 
high-resolution imager (HRI) reveals an ISM pressure of 
px ~ 5 x 10~ 10 dyncm -2 (C. B. Peres, private communica- 
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tion). Thus the highest minimum pressures obtained from 
radio observations exceed the ISM pressure by a factor of 2. 
The amount by which the radio halo is overpressured with 
respect to the ISM can be estimated by considering the ex- 
pansion of an overpressured bubble. Bicknell & Begelman 
(1996) present such a model for the M87 system and show 
that the inner radio halo should be overpressured with re- 
spect to the ISM by a factor of ~ 3. Independent support for 
this mild overpressure comes from the observed optical emis- 
sion line filaments discovered by Ford & Butcher (1979) and 
studied in greater detail by Sparks, Ford and Kinney (1993). 
Bicknell & Begelman (1996) show how the expansion of this 
mildly overpressured halo energizes the optical emission line 
filaments. 

To conclude this argument, the pressure of the inner 
radio halo is likely to be p w 1.5 x 10 -9 dyncm -2 , leading 
to a stored energy within the halo of 6 x 10 56 erg. The total 
kinetic luminosity required to supply this energy in 10 6 yr 
is 2 x 10 43 ergs _1 . If there is a symmetric counter jet (as 
assumed in the classic twin-jet model) , the single-jet kinetic 
luminosity which contributes to the relativistic particle en- 
ergy in the halo is 1 x 10 43 ergs -1 . 

4-1.3 Time-averaged kinetic luminosity, Lk 

Summing the above contributions, the kinetic luminosity of 
the jet (averaged over the duration of the current period of 
activity) is Lk = 10 43 Z/43 erg s _1 where L43 ~ 1. 

An independent argument against a high-power jet can 
be put forward by considering knot- A. In the kinematic 
model of Biretta (1993), knot-A represents a strong trans- 
verse shock in the flow where the jet decelerates from r ~ 3 
to T ~ 1.5. A substantial part of the kinetic luminosity 
would then be converted to other forms at the shock. Thus, 
the fact that we only see ~ 2 x 10 42 ergs -1 radiated at 
this knot suggests that the jet does not possess a kinetic 
luminosity that is orders of magnitude above this value. 
Note that the fact that the flow remains collimated after 
knot-A suggests that the kinetic energy is not transformed 
purely into internal energy of the jet material. Even in the 
more physically plausible oblique shock models of Bicknell 
& Begelman (1996), a jet with a high kinetic luminosity 
(Lk ^ 10 44 erg s _1 ) would be expected to radiate more than 
is observed. Thus, a high jet power necessarily entails prob- 
lems in understanding the properties of knot-A. 

Of course, the above arguments cannot rule out the 
possibility that the jet has undergone a recent dramatic in- 
crease in kinetic luminosity. Such an increase could lead to 
a powerful VLBI core whilst still keeping the total (aver- 
age) energetics of the system at the level observed. If the 
jet has undergone a dramatic increase in power, one might 
expect to observe the effects of this within the inner jet (i.e. 
between the core and knot-A, where we have independent 
arguments for a Lk ~ 10 43 ergs -1 jet.) However, the inner 
jet is fairly featureless, displaying only the relatively weak 
knots D, E, and F. Furthermore, there is no precedent for the 
core of M87 to undergo large amplitude (i.e. order of magni- 
tude) variability. Thus, to postulate a transient increase in 
jet power would be an ad-hoc way to bypass the energetic 
constraints presented below. 



4.2 Lk and the physical jet parameters 

To relate the limits on Lk to physical parameters of the jet, 
we assume the jet to be uniform in the sense that the density 
and velocity are constant across a cross-section of the jet at 
a given distance, Z, along the jet. If positive charge carriers 
have mass m+ per unit electron charge, the total jet kinetic 
luminosity is given by 

L K = nr(Z) 2 (3(r - l)rn(m c ( 7 ) + m+( 7+ ))c 3 , (10) 

where r(Z) is the physical radius of the jet a distance Z 
away from its origin. The inclusion of the (7) and (7+) terms 
allow for energy that is advected within the bulk flow for the 
electron and positive charge carrier populations respectively. 

From the above equation, n can be expressed in terms 
of quantities for which we have observational constraints 
and evaluated for the case of an electron-proton jet (m+ = 
1836m c ) and an electron-positron jet (m+ = m c ) at the 
5 GHz self-absorption point. We assume that L43 = 1, T = 3 
and (7) = 10 (corresponding to the 7 m in ~ 1 case). Relax- 
ation of these assumptions are discussed in the next section. 
For the electron-positron jet (m+ = m e , (7+) = (7) ^ 10) 
the result is n ~ 1.3 x 10 2 cm -3 . For the electron-proton 
jet, the proton population will be essentially cold unless 
the average energy per electron is comparable to or greater 
than the proton rest mass (i.e. (7) > 2000). In this case 
(m+ = 1836ra c , (7+) 1) we deduce n « 1.4cm -3 . These 
results are shown as horizontal lines on Fig. 1. 

Examination of Fig. 1 strongly suggests the jet to 
be dominated by an electron-positron plasma rather than 
an electron-proton plasma. The 22 GHz VLBI data and 
100 GHz VLBI data lead to a similar conclusion. This is 
the main result of the present paper. 

The only way for the above constraints to be compatible 
with an electron-proton jet is if the relativistic electron pop- 
ulation has a low-energy cutoff at jminin c c 2 with 7> 1. Us- 
ing either the 5 GHz or 22 GHz constraints requires 7 m i n > 
10. The 100 GHz constraints require 7 m i n > 100 (assum- 
ing the core is still self-absorbed at this frequency). Phys- 
ically, increasing 7 m i n allows the normalization of the rela- 
tivistic particle distribution, No, to increase whilst holding 
the density, n, fixed. Since the synchrotron self-absorption 
constraints depend only on No , varying 7 m i n allows the syn- 
chrotron constraints to decouple from the kinetic energy con- 
straints. In particular, line B on Fig. 1 will move downwards 
as 7 m in increases such as to keep the combination n7 m i n B 2 
constant. Observational and theoretical constraints on 7 m i n 
will be considered in the next Section. It will be seen that the 
range of parameter space available for an electron-positron 
jet is relatively small and not easily justifiable on physical 
grounds. Furthermore, the electron-proton solution is vul- 
nerable to any future tightening of these observational con- 
straints. 



5 DISCUSSION 

The above constraints argue for an electron-positron jet as 
opposed to an electron-proton jet. In this section we examine 
the robustness and self-consistency of this result. We then 
discuss some astrophysical implications. 
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5.1 The dimension of the self-absorbed core 

We have taken the angular size of the resolved core, as de- 
termined from 5 GHz VLBI, to represent the diameter of the 
jet at the point where it just becomes self-absorbed. More 
precisely, this should be considered an upper limit to the jet 
dimension. The VLBI core dimension might plausibly rep- 
resent the (projected) length of the entire self-absorbed jet. 
The effect of this uncertainty on our constraints is simply 
addressed by re- calculating Fig. 1 for a smaller jet diameter. 
Suppose the physical jet diameter was a factor of 2 smaller 
than that inferred from VLBI observations. The synchrotron 
surface brightness constraint, eqn. (6), becomes tighter (i.e. 
the upper limit on B reduces) by a factor of 4 (since the 
surface brightness is twice that observed). The optical depth 
constraint also become tighter by a factor of 2. The net effect 
is to raise the (synchrotron) lower limit on n7 m i n by a factor 
of 32. By constrast, the values of n derived from the kinetic 
luminosity will only be raised by a factor of 4. Thus, the net 
effect is to tighten the argument against an electron-proton 
jet. 

If the jet in the core is, in fact, significantly narrower 
than the VLBI resolution then even the electron-positron 
solution violates the synchrotron constraints. This is an ar- 
gument against such a narrow jet. 



5.2 The effect of a high velocity jet-core 

It is physically plausible that the jet starts as a high Lorentz 
factor beam (r 3> 10) and subsequently decelerates via dis- 
sipative processes or entrainment of surrounding material 
(Begelman, Rees & Sikora 1994; Sikora, Begelman & Rees 
1994). The VLBI jet may then possess a high Lorentz factor 
core surrounded by a slower sheath (F ~ few). The emis- 
sion from the core would be largely beamed out of our line 
of sight and thus would be swamped by emission from the 
slower sheath. Our above synchrotron constraints would ap- 
ply only to the outer sheath since it would dominate the 
observed radio emission. However, since the kinetic energy 
constraints are based on the large scale energetics of the 
system, postulating the presence of an unseen rapid core re- 
duces the inferred kinetic luminosity of the outer sheath and 
thus makes the case against an electron-proton jet stronger. 



5.3 Synchrotron self-Compton radiation and 
observational limits on 7 m i n 

5.3.1 Synchrotron self-Compton radiation 

The production of high-energy synchrotron self-Compton 
(SSC) radiation is an unavoidable consequence of the syn- 
chrotron process. We must check that our jet models are 
consistent with the constraint that the SSC radiation does 
not violate the observed X-ray limits on the core flux. We 
can also use this constraint to place observational limits on 

7min • 

The SSC X-ray flux density, Fx from the self-absorbed 
region (size r) is given by G92 as 



Fx = 



2aF r 



t(a) 



V Vx ) \ V m / 



in the synchrotron spectrum and t(a) is a function tabulated 
in G92. We typically assume that Vh lies somewhere between 
optical and X-ray frequencies. Evaluating this flux density 
at 1 keV for parameters relevant to the 5 GHz self absorbed 
region gives 



F x « 3.0 x 10 n n7 mi „ In {v h /v m ) Jy. 
where 

ln(j^/z^ m ) ~ 10 



(12) 



(13) 



(11) 



where vx is the X-ray frequency, v\> is the high-energy cutoff 



for physically plausible values of v\,. For consistency, this 
must be less than the observed X-ray core flux density of 
~ 3.5 x 10~ 7 Jy (BSH). We note that using the result of 
Section 3.3 we can rewrite this as a function of the VLBI 
core flux density and the magnetic field B only. This gives 
a lower limit on the magnetic field of B > 0.01 G. 

For pair jets, we have shown from our kinetic lumi- 
nosity arguments that n w 10 2 cm' 3 . The production of 
synchrotron self-Compton radiation is then consistent with 
observations provided 7 m i n < 10 (in which case, a substan- 
tial fraction of the observed X-rays will originate from the 
SSC process). For an electron-proton jet (n fa 1cm" 3 ) the 
corresponding limit is 7 m i n < 10 3 . 



5. 3. 2 SSC limits and 7mm 

The SSC limits give B > 0.01 G. Thus, an electron en- 
ergy distribution that cuts off at 7 m i n will produce a low- 
energy cutoff in the synchrotron spectrum at frequency 
u c ~ 7 min (B/G) MHz. Observationally, u c < 1 GHz leading 
to the condition that 7 m i n < 100. Thus, if 10 < 7 m in < 100, 
the electron-proton jet is consistent with our 5 GHz and 
22 GHz constraints. If the core emission is still dominated 
by self-absorbed radiation at 100 GHz, the corresponding 
VLBI constraints are incompatible with an electron-proton 
jet since compatibility would demand 7 m i n ~ 100 which 
would lead to an unobserved low-frequency cut off at v c ~ 
100 GHz (note that the SSC limit on the magnetic field at 
the 100 GHz self-absorption location is B > 10 G) . 



5.4 Theoretical issues related to 7 m i n 

In the absence of further observational limits on 7 m in, we 
briefly discuss some related theoretical issues. 

Let us suppose that the relativistic electron population 
at the 5 GHz self-absorbed point did, in fact, possess a low- 
energy cutoff corresponding to 7 m i n ~ 10-100. Such an en- 
ergy distribution cannot be a relic of the conditions imposed 
near the base of the jet: synchrotron and inverse Compton 
losses would have cooled a significant number of these elec- 
tron to energies below 10m c c 2 . Thus, some in-situ physical 
mechanism must be maintaining this low-energy cutoff (i.e. 
by either reheating the electrons/positrons or injecting fresh 
high-energy pairs). 

One mechanism that has been proposed for injecting 
high-energy pairs into a jet is that of high-energy hadronic 
(primarily proton-proton, hereafter p-p) collisions (Falcke & 
Biermann 1995; Biermann, Strom & Falcke 1995). Such col- 
lisions produce both charged and neutral pions. The charged 
pions then decay into electrons or positrons, corresponding 
to the n~ and tt + cases respectively, with an initially high 
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energy (7 ~ 100) due to the large rest mass ol the pions. 
However, the p-p interaction cross section is approximately 
ffpp w (jt/10 (and decreases as 7 increases). Thus, the opti- 
cal depth to p-p collisions at the location where our VLBI 
constraints are applied is < 10~ 8 (where we have used the 
fact that the number density of protons is constrained to 
be n p < 1cm" 3 by the kinetic luminosity constraint). We 
conclude that this mechanism is not relevant in this case. 
In fact, it would only be relevant at the base of very pow- 
erful (and compact) jets. We note that hadronic processes 
may be of importance in compact plasmas near the centre 
of the accretion flow (Sikora et al. 1987; Begelman, Rudak 
& Sikora 1990). In the absence of any other mechanism that 
inject high-energy electrons, any low-energy cutoff seen at 
the 5 GHz self-absorption point is not likely to be due to 
such injection. 

Various heating mechanisms that could lead to 7 m i n > 1 
have been discussed by G92. The most likely heating mecha- 
nism is that due to synchrotron self-absorption (the so-called 
'synchrotron boiler'; Ghisellini, Guilbert & Svensson 1988). 
For the region under consideration here, electrons radiating 
at the self-absorpion frequency (5 GHz) have 7 ~ 100. If 
only synchrotron processes were important (i.e. the plasma 
is magnetically dominated), the electron energy distribution 
below this would possess a quasi-thermal component (due to 
the efficient excahnge of energy between the relativistic elec- 
trons). The effective low-energy cutoff would then lie some- 
where below 7min ~ 100. However, if inverse Compton losses 
are important (as is probably the case for jets on the VLBI 
scale), the effective low-energy cutoff will be much reduced 
and, indeed, synchrotron heating can become irrelevant. A 
detailed treatment of this phenomenon is beyond the scope 
of this paper. 

5.5 Annihilation of the pairs 

The cross-section for annihilation of cold (i.e. non- 
relativistic) pairs is given by cr a = (3/8)<7t where or is 
the Thomson cross-section. [As the pairs become relativis- 
tic, the cross-section declines in a manner analogous to the 
Klein-Nishina decline of the Thomson cross-section.] In the 
comoving frame, the rate of annihilation of cold pairs is then 

h — -ccttw , (14) 
8 

leading to an annihilation timescale of 
8 

tann — ~ • (15) 

ScaTn 

In the outer, low-density parts of the jet this timescale is long 
and annihilation is not an important process. However, in 
the higher densities found near the jet core the annihilation 
rate increases. It is instructive to estimate the distance from 
the nucleus at which the annihilation timescale, t ann equals 
the flow timescale tfl ow = Z/Tc where Z is the distance from 
the core. If the jet has a constant opening angle and constant 
velocity so that the density falls off as n oc Z~ 2 , equating 
tann and tfl ow gives the annihilation radius 

^ann — — , (.10) 

where n is the particle density at a distance Z along 
the jet. Our (cold-jet) kinetic luminosity constraints give 



nr = 4 x 10 cm" . If we assume a half opening angle of 
> 0.1 rad (such as to agree with the opening angle on larger 
scales) then Z < lOr. Together with V = 3, this evaluates 
to give Z = 3 < 10 13 cm. This is to be compared with the 
Schwarzschild radius of R s « 1 x 10 15 cm (M = 3 x 10 9 M 
assumed) . Thus it is feasible for the jet to have originated 
from a pair plasma formed at the centre of the accretion 
flow. 

5.6 The FR-I/FR-II dichotomy 

Celotti & Fabian (1993) utilized SSC constraints on a sam- 
ple of radio galaxies and radio-loud quasars (dominated 
by FR-II type objects) and argued for either an electron- 
positron jet with 7 m in ~ 1 or an electron-proton jet with a 
7min ~ 100. For these powerful sources, it was found that 
the annihilation radius lies well outside the region occupied 
by any compact pair plasma associated with the inner accre- 
tion flow. Thus it is not possible to have the pair jet flowing 
freely from the immediate vicinity of the black hole. This 
argument was used to favour the electron-proton case. How- 
ever, we have presented arguments for an electron-positron 
jet in M87 (a classical FR-I source). These results suggest 
one of two possibilities: 

a) FR-I and FR-II may possess jets of an intrinsically 
different nature, at least by the time the jet has propa- 
gated to VLBI scales: i.e. FR-I sources may possess electron- 
positron VLBI jets whereas FR-II sources may possess VLBI 
electron-proton jets. The morphological differences between 
FR-I/FR-II sources may reflect this underlying difference in 
the jet (for example, via the ease of collimation) . However, it 
is difficult to envisage the physical processes leading to such 
a dichotomy and how they could be related to the source 
power. 

b) The annihilation radius constraint does not apply 
and both classes of sources possess electron-positron jets 
with 7 m in ~ 1. In this case, the jet would have to propa- 
gate from the vicinity black hole out to at least the anni- 
hilation radius in some other form such as a Poynting flux 
dominated flow. Beyond the annihilation radius, electromag- 
netic cascades can lead to copious pair production within the 
Poynting flux dominated flow. 



6 CONCLUSIONS AND DEFINITIVE 
OBSERVATIONS 

The M87 jet has the best determined geometry and kine- 
matics of any extragalactic jet in the sky. M87 is also one 
of the very few galaxies in which the mass of the central 
black hole (2.4 ± 0.7 x 10 9 M ) is well constrained (by HST 
observations). This makes M87 an obvious system in which 
to study the physical properties of extragalactic jets. 

We have utilized the standard theory of synchrotron 
self-absorbed cores in order to constrain the magnetic field, 
B, and the relativistic (proper) particle density, n, of the 
jet. The 5 GHz data implies n7 m i n ^ 10 cm" 3 . For this jet to 
carry the kinetic liminosity inferred from global energetic ar- 
guments, the density is either n ~ 1 cm -3 (electron-proton 
jet) or n ~ 10 2 cm" 3 (electron-positron jet). We cite this 
as evidence for an electron-positron dominated jet rather 
than an electron-proton dominated jet. An electron-proton 
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jet with 10 < 7 m in < 100 is consistent with the present (con- 
servative) constraints but is extremely vulnerable to further 
tightening of the VLBI limits. If the core seen in 100 GHz 
VLBI data is self-absorbed, the electron-proton jet is ruled 
out. Some theoretical arguments are also presented detailing 
the problems of the high-7 m i n scenario. 

To make further observational progress, the self- 
absorbed core has to be identified and resolved. This re- 
quires high-resolution multi-frequency VLBI observations. 
The sub-arcsec X-ray imaging capability of AXAF will sep- 
arate the core X-ray emission from other components (e.g. 
inner knot emission) and thus allow stricter limits to be 
placed on the SSC radiation of the self-absorbed core. Such 
data, together with arguments similar to those presented 
here, should allow the matter content of this archetypal ex- 
tragalactic jet to be unambiguously determined. 
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